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A multi-cluster model of traffic flow is studied, in which the motion of cars is described by a 
stochastic master equation. Assuming that the escape rate from a cluster depends only on the 
cluster size, the dynamics of the model is directly mapped to the mathematically well-studied zero- 
range process. Knowledge of the asymptotic behaviour of the transition rates for large clusters 
allows us to apply an established criterion for phase separation in one-dimensional driven systems. 
The distribution over cluster sizes in our zero-range model is given by a one-step master equation in 
one dimension. It provides an approximate mean-field dynamics, which, however, leads to the exact 
stationary state. Based on this equation, we have calculated the critical density at which phase 
separation takes place. We have shown that within a certain range of densities above the critical 
value a metastable homogeneous state exists before coarsening sets in. Within this approach we have 
estimated the critical cluster size and the mean nucleation time for a condensate in a large system. 
The metastablity in the zero-range process is reflected in a metastable branch of the fundamental 
flux-density diagram of traffic flow. Our work thus provides a possible analytical description of 
traffic jam formation as well as important insight into condensation in the zero-range process. 

PACS numbers: 89.40.-a, 02.50.Ey, 64.75+g 



I. INTRODUCTION 

The development of traffic jams in vehicular flow is 
an everyday example of the occurence of phase separa- 
tion in low-dimensional driven systems, a topic which 
has attracted much recent interest (see, e. g., 0,0 and 
references therein). In Q the existence of phase sepa- 
ration is related to the size-dependence of domain cur- 
rents and a quantitative criterion is obtained by consid- 
ering the zero-range process (ZRP) as a generic model 
for domain dynamics. Phase separation corresponds to 
the phenomenon of condensation in the ZRP (see Q for 
a recent review) whereby a macroscopic proportion of 
particles accumulate on a single site. 

In this paper we use such a zero-range picture to study 
the phase separation in traffic flow. A stochastic master 
equation approach in the spirit of the probabilistic de- 
scription of transportation [5( enables us to calculate the 
critical parameters of our model. We pay particular at- 
tention to the initial nucleation of a condensate in the 
ZRP, complementing previous work 0, Q on the late- 
stage cluster coarsening. Significantly, we find that prior 
to condensation the system can exist in a homogeneous 
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metastable state and we provide estimates of the criti- 
cal cluster size and mean nucleation time. Finally, we 
apply these results to the description of traffic, obtain- 
ing a fundamental flux-density diagram which includes 
a metastable branch. Metastability and hysteresis ef- 
fects have been observed in real traffic, see, e. g., d, 
for discussion of empirical data and the various different 
modelling approaches. For previous work focusing on 
the description of jam formation as a nucleation process, 
see [HA 111) . 

II. THE MODEL 

We consider a model of traffic flow, where cars are 
moving along a circular road. Each car occupies a cer- 
tain length of road I. We divide the whole road of total 
length L in to cells of size £. Each cell can be either empty 
or occupied by a car, just as in cellular automaton traf- 
fic models (see, e. g., [3, 13 an d references therein). Most 
such models use a discrete-time update rule, for example, 
see 01 f° r a class of traffic models related to a parallel 
updating version of the ZRP. In contrast, we consider 
the development of our system in continuous time. The 
probability per unit time for each car to move to the 
next cell is given by a certain transition rate, which de- 
pends on the actual configuration of empty and occupied 
cells. This configuration is characterised by the cluster 
distribution. An uninterrupted string of n occupied cells, 
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bounded by two empty cells, is called a cluster of size 
ri. The clusters of size n = 1 are associated with freely 
moving cars. The first car in each cluster is allowed to 
move forward by one cell. The transition rate w n of this 
stochastic event depends on the size n of the cluster to 
which the car belongs. In this case w\ is the mean of the 
inverse time necessary for a free car to move forward by 
one cell. The transition rate w\ is related to the distribu- 
tion of velocities in the free flow regime or phase, which 
is characterised by a certain car density Cf reo . For small 
densities, expected in the free flow phase in real traffic, 
the interaction between cars is weak and, therefore, the 
transition rate w± depends only weakly on the density 
Cfree- Hence, in the first approximation we may assume 
that wi is a constant. 

This model can be directly mapped to the zero-range 
process (ZRP). Each vacancy (empty cell) in the origi- 
nal model is related to a box in the zero-range model. 
The number of boxes is fixed, and each box can contain 
an arbitrary number of particles (cars), which is equal 
to the size of the cluster located to the left (if cars are 
moving to the right) of the corresponding vacancy in the 
original model. If this vacancy has another vacancy to 
the left, then it means that the box is empty. Since the 
boundary conditions are periodic in the original model, 
they remain periodic also in the zero-range model. In 
this representation, one particle moves from a given box 
to the right with transition rate w n , which depends only 
on the number of particles n in this box. In the grand 
canonical ensemble, where the total number of particles is 
allowed to fluctuate, the stationary distribution over the 
cluster-size configurations is the product of independent 
distributions for individual boxes. The probability that 
there are just n particles in a box in a homogeneous phase 
is E1II1 P(n) oc z n /Y[ n m=l w m for n > 0, P(0) being 
given by the normalisation condition. Here z = e fl / kBT 
is the fugacity - a parameter which controls the mean 
number of particles in the system. 
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III. MASTER EQUATION 

This result can be obtained and interpreted within the 
stochastic master equation approach |4|. Assuming the 
statistical independence of the distributions in different 
boxes, we have a multiplicative ansatz 



dP(n, t) 
di 

dP(0,t) 



dt 



(w)P(n - 1, t) + w n+1 P(n + 1, t) 
[(w)+w n ]P(n,t) : n>l, (2) 
Wl P(l,t) - (w)P(0,t) , (3) 



where 



(w)(t)=^2w k P(k,t) 



(4) 



fc=i 



is the mean inflow rate in a box. The ansatz is an 
exact property of the stationary state of the grand canon- 
ical ensemble or, alternatively, of an infinitely large sys- 
tem ^3|. Hence, in these cases, the master equations © 
and @ give the exact stationary state while providing a 
mean- field approximation to the dynamics of reaching it. 

The stationary solution P{n) corresponding to 
dP(n,t)/dt = can be found recursively, starting from 
n = 0. It yields the known result [HT^ . il ^ | 



- 1 

p(n)=p(p)(w)» n — 



(5) 



for n > 0, where P(0) is found from the normalisation 
condition. 

Denoting by M the number of boxes, which corre- 
sponds to the number of vacancies in the original model, 
the mean number of cars on the road is given by (iV) = 
M(rt), where 



(n> 



nP(n) 



(6) 



is the average number of particles in a box. Note that in 
the grand canonical ensemble the total number of cars as 
well as the length of the road L fluctuate. For the mean 
value, measured in units of £, we have (L) = M + (N). 
Hence, the average density of cars is 



(TV) 
(L) 



(n> 



l + (n> 



(7) 



According to 0, @, and JSJ, we have the following 
relation 



P 2 (k,m,t) = P{k,t) P{m,t) 



(1) 



for the joint probability P2{k, to, t) that there are to par- 
ticles in one box and k particles in the neighbouring 
box on the left at time t. This approximation leads to 
the mean-field dynamics described by the master equa- 



e n( w ) n n — 

c _ 7t=l m— 1 

i+EH n n — 

n—1 m—1 u^rn 



(8) 



from which the stationary mean inflow rate (w) can be 
calculated at a given average density c. 
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FIG. 1: {w)/woo vs density c at b — 1 for different a. 



IV. TRANSITION RATES AND PHASE 
SEPARATION 

Now we make the following choice for the transition 
rate dependence on the cluster size n: 



1 



for n > 2 



(9) 



the value of w± being given separately, since it is re- 
lated to the motion of uncongested cars, whereas w n 
with n > 2 represents the escaping from a jam of size 
n. Although an individual driver does not know how 
many cars are jammed behind him, the effective current 
of cars from a jam, represented by w n , is a collective ef- 
fect which is expected to depend on the correlations and 
internal struture (e. g., distribution of headways) within 
the cluster 0. A monotonously decreasing dependence 
on cluster size, such as (0, can be considered as a type 
of slow-to-start rule — the longer a car has been station- 
ary the lar ger the probability of a delay when starting 
(cf. |ll EH 111 E]|). 

We now explore the consequences of the choice (J5J in 
terms of the ZRP phase behaviour and its implications 
for the description of traffic flow. In numerical calcula- 
tions we have assumed Woo = 1/roo = 1 and w\ = 5, by 
choosing the time constant T m as a time unit, whereas 
the control parameters 6 and a have been varied. 

If cr > 1, as well as for b < 2 at a = 1, Eq. © has a 
solution for any density < c < 1 (see dashed and dot- 
ted curves in Fig.^l. This implies that the homogeneous 
phase is stable in the whole range of densities, i. e., there 
is no phase transition in a strict sense. If cr < 1 (solid 
curve in Fig. [IJ, as well as for 6 > 2 at a = 1, (w)/woo 
reaches 1 at a critical density < c cr < 1 , and there is no 
physical solution of iJSJ for c > c cr . This means that the 
homogeneous phase cannot accommodate a larger den- 
sity of particles and condensation takes place at c > c cr . 

This behaviour underlies the known criterion for phase 
separation in one-dimensional driven systems [3j . For il- 
lustration, we comment that in the multi-cluster model 
considered in [19| the transition rates do not depend on 



the cluster sizes, only the inflow rate in a cluster depends 
on the overall car density and fraction of congested cars. 
This corresponds to the case 6=0, where, according 
to the criterion, no macroscopic phase separation takes 
place - in agreement with the theoretical conclusions and 
simulation results of [3- In contrast, a class of mi- 
croscopic models was introduced in I20l | where correla- 
tions within the domain (jam) give rise to currents of the 
form with a = 1 and b > 2; phase separation is then 
observed. 

At c < c cr m our model the cluster distribution func- 
tion P(n) decays exponentially fast for large n whereas 
the decay is slower at c = c cr . It is well known that 
the decay in this case is power-like for a = 1, i. e., 
P{n) ~ n~ b 0. For < a < 1, we find that 
the leading behaviour is stretched exponential, i. e., 
P(n) ~ exp [— 6n 1_cr /(l — cr)] , in agreement with the 
result stated in j4|. Within our approach we have also 
derived the sub-leading terms which turn out to be rel- 
evant for < cr < 1/2. This calculation is presented in 
the Appendix and further illustrates the rich behaviour 
of the model as cr is varied. 

At (w) I Woo = 1 the inflow (w) in a macroscopic cluster 
of size n — ► oo is balanced with the outflow u>oo. This 
means that at overall density c > c cr the homogeneous 
phase with density c cr is in equilibrium with a macro- 
scopic cluster, represented by one of the boxes containing 
a non- vanishing fraction of all particles in the thermo- 
dynamic limit [a, 0|- Hence, (w)/w ao — 1 holds in the 
phase coexistence regime at c > c cr . 

According to JJJ), the critical density c cr is given by 



(10) 



1 + <n)c 



where (ri) cr is the mean cluster size at the critical density. 
Since (w) = holds in this case, we have 



(n)c 



oo n ^ 

e nw ^o n — 

n— 1 m=l ^rn 

oo n 1 

i+e^ n — 

n=l m=l w m 



(11) 



The critical density, calculated numerically from l|10|) 
and I jllj l as a function of parameters a and b, is shown in 
Figs. |21 and |3 respectively. In contrast to the situations 
discussed previously in the literature, in our model wi 
is not given by the general formula @ but is an inde- 
pendent parameter. This distinction leads to quantita- 
tively different results, for example, the critical density 
for cr = 1 is analytically shown to be [22| 



6(6+1) 



(6-1)[2(6 + 1) + ib 1 (6-2)] 



(12) 
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FIG. 2: Critical density as a function of control parameter a 
for different values of b. 




FIG. 3: Critical density as a function of control parameter b 
for different values of a. 



V. METASTABILITY 

Suppose that at the initial time moment t = the 
system is in a homogeneous state with overall density 
slightly larger than c cr . Here we study the development 
of such a state in the mean-field approximation provided 
by (J2J and 0. With this initial condition, the mean 
inflow rate in a box (w) is slightly larger than that at c = 
c cr , i. e., (w) = Woo + £ holds with small and positive s. 
Hence, only large clusters with w n < Woo+e have a stable 
tendency to grow, whereas any smaller cluster typically 
(except a rare case) fluctuates until it finally dissolves. 
In other words, the initially homogeneous system with no 
large clusters can stay in this metastable supersaturated 
state for a long time until a large stable cluster appears 
due to a rare fluctuation. 

Neglecting the fluctuations, the time development of 
the size n of a cluster is described by the deterministic 
equation 



dn 
~dt 



(w) 



w„ 



(13) 



with n < n cr tend to dissolve, whereas the overcritical 
ones with n > n cr tend to grow, where the critical cluster 
size n cr is given by the condition 



(w) = w nc 



Using J3J yields 



[w)/w 



1 



l/o 



(14) 



(15) 



In this rounded to an integer value. 

This deterministic approach describes only the most 
probable scenario for an arbitrarily choosen cluster of a 
given size. It does not allow one to obtain the distribu- 
tion over cluster sizes: the deterministic equation 113|) 
suggests that all clusters shrink to zero size if they are 
smaller than n cr at the beginning, whereas the real size 
distribution arises from the competition between oppo- 
site stochastic events of shrinking and growing. Assum- 
ing that the distribution of relatively small clusters con- 
tributing to (n) is quasi-stationary, i. e., that the detailed 
balance (equality of the terms in J5J and 10 describing 
opposite stochastic events) for these clusters is almost 
reached before any cluster with n > n cr has appeared, 
we have 



i)~£nP(n) 



(16) 



for such a metastable state. In this case from (JJJ) we 
obtain 



n cr n 1 

£ n (w) n n — 

£ n—1 m—1 ^rn 



1 _ c ncr n I 

i+EW" n — 

n— 1 m—1 ^>m 



(17) 



According to this equation, the undercritical clusters 



instead of JHJ for calculation of (w) in this homogeneous 
metastable state. The critical cluster size is found self 
consistently solving (|15fl and (|17fl as a system of equa- 
tions. From l|15|l we see that the critical cluster size n cr 
diverges at c — > c cr , since (w) — * w^,. The results of cal- 
culation of n cr (with rounding down to an integer value) 
at a = 0.5, 6=1 and at a = 1, b — 3 are shown in Figs. 01 
and [S] respectively. 



VI. THE MEAN NUCLEATION TIME 

To our knowledge, the nucleation time, i. e., the mean 
time the system stays in the homogeneous metastable 
state, has not been considered up to now for the zero- 
range process. A characteristic time scale for a particle 
to escape from a jam has been evaluated in for the 
related bus route model. This time scale, however, refers 
to already formed finite-size jams and is therefore differ- 
ent to the nucleation time we are interested in. 
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c cr 0.6 



FIG. 4: Critical cluster size vs density at a = 0.5 and 6=1. 
The critical density is indicated by a vertical dashed line. 



given by 

M-l 

V M (t) = MV(t) x 



Z 

1 - j V{t')dt' 
o 



MPihcx]> j -M J V(t')dt' j . ( IN) 

The latter equality holds for large M, since all M boxes 

t 

are equivalent, and therefore the probability JP(t')dt' 

o 

that the nucleation occurs in a given box within a char- 
acteristic time interval t ~ (T)m is a small quantity of 
order 1/M. The mean nucleation time for the system of 
M boxes is 



lOOOr 

n 

cr 
500- 



L 




cr 



0.5 



0.6 



FIG. 5: Critical cluster size vs density at cr = 1 and 6 = 3. 
The critical density is indicated by a vertical dashed line. 



Within the framework of mean-field dynamics, the 
mean nucleation time in our model can be evaluated 
as follows. Let "P(i) be the probability density of the 
first passage time to exceed the critical number of par- 
ticles n cr in a single box. By our definition, the nucle- 
ation occurs when one of the M boxes reaches the clus- 
ter size n cr + 1. The probability that it occurs first in a 
given box within a small time interval [i; t + dt] is thus 
t 



V{t)dt x 



1 



JV(t')df 
o 



M-l 



according to our assump- 



The 



tion that the boxes are statistically independent. 

r t -iM-l 
term 1 — J V{t')dt' is the probability that in all 

o J 

other boxes, except the given one, the overcritical clus- 
ter size n cr + 1 has still not been reached. Since, the 
nucleation can occur in any of M boxes, the nucleation 
probability density Vuify for the system of M boxes is 



(t) m = J tr M (t)dt . 

o 



(19) 



Here (T) i is the mean first passage time for a single box. 

In order to estimate (T)m according to (|T%)l and ltH?|l . 
one needs some idea about the first passage time proba- 
bility density for one box V(t). This is actually the prob- 
lem of a particle escaping from a potential well. Since we 
start with an almost homogeneous state of the system, 
we may assume zero cluster size n = as the initial con- 
dition. The first passage time probability density can be 
calculated as the probability per unit time to reach the 
state n cr + l, assuming that the particle is absorbed there. 
It is reasonable to assume that after a certain equilibra- 
tion time t eq , when a quasi-stationary distribution of the 
cluster sizes within n < n cr is reached, the escaping from 
this region is characterised by a certain transition rate 
w esc . Hence, for t > t eq we have 



V(t) 



1 - J V{t')dt' 



(20) 



where the expression in square brackets is the probability 
that the absorption at n cr + 1 has still not occured up 
to the time moment t. At high enough potential barriers 
(large mean first passage times) the short-time contribu- 
tion to the integral is irrelevant and, by means of i|19|l . 
the solution of l|2U|l can be written as 



where 



T(t) 



1 



(T) 



■ exp 



(21) 



Obviously, this approximate solu- 



tion of the first passage problem is not valid for very 
short times t <^it eq , since the short-time solution should 
explicitly depend on the initial condition. In particu- 
lar, if we start at n = 0, then the state n cr + 1 cannot be 
reached immediately, so that 'P(O) = 0. Nevertheless Q21JI 
can be used to estimate the mean nucleation time {T)m 
provided that (T)m > teq- 
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5x10 



1 













eq 



FIG. 6: Comparison between the theoretical approxima- 
tion l|21|l for V(t) (smooth curves) and mean-field simula- 
tion results (fluctuating curves) shown in a longer (top) and 
in a shorter (bottom) time scale. The vertical dashed line 
indicates the equilibration time t eq ~ 500 for the set of pa- 
rameters (7 = 0.5, 6 = 1, and c = 0.84 represented by the 
upper curves in both pictures. The other curves correspond 
to a = 1, b = 3, and c = 0.61. 



FIG. 7: Mean nucleation time vs density at a — 0.5, 6 = 1, 
and M = 10 6 . The critical density is indicated by a vertical 
dashed line. 



By inserting i|21|) in i|19fl we obtain 



(T) 



M 



M 



OO 

(T)i J x e~ x exp (-M [l - e~ x ] ) dx (22) 



after changing the integration variable tj {T)i — > x. Tak- 
ing into account that only the region x ~ 1/M con- 
tributes to the integral at large M, we arrive at a very 
simple expression 



(T) 



M 



M 



(23) 



relating the mean first passage time or nucleation time 
in a system of M boxes with that of one box. The latter 
can be calculated easily by the known formula [23| 



(24) 



n=0 



m—0 



We have checked the correctness of these theoretical 
expectations within the mean-field dynamics represented 
by and J31 by comparing them with the results of sim- 
ulation of stochastic trajectories generated according to 
these equations. The simulation curves for V(t) at two 
different sets of parameters: a = 0.5, b = 1, c = 0.84 
(with n cr — 48), and a = 1, b = 3, c = 0.61 (with 
n cr = 35) are shown in Fig.Elin two different time scales. 
As we see, Eq. (|21H is a good approximation for large 
enough times t > t eq . For definiteness, we have iden- 
tified the equilibration time t eq with the crossing point 
of the theoretical and simulated curves. An interesting 
additional feature is the presence of an apparent nucle- 
ation time lag, which is about ti ag « 60 for the first set 
of parameters and about ti ag ~ 30 for the second one. 
Evidently, the first passage time probability density V(t) 
tends to zero very rapidly when t decreases below ti ag . 



where P(0) = 1 and P(n) = ]J ((w)/wk) with n > 1 

represent the unnormalised stationary probability distri- 
bution. 

The mean nucleation time versus density c, calculated 
from (ESI and |22J) at M — 10 6 is shown in Figs. □ and El 
Fig. [prefers to the case a = 0.5 and b = 1, whereas Fig- El 
is for the case a = 1 and 6 = 3. These figures show that 
the mean nucleation time increases dramatically as the 
critical cluster size n cr increases (see the corresponding 
plots in Figs. 01 and El approaching the critical density 
c cr . 

According to our previous discussion, estimate (|23|l is 
valid for large enough mean nucleation times {T)m > 
t eq , in particular, when approaching the critical density 
c \ c cr at any large but fixed M. It is not valid in 
the thermodynamic limit M — > oo at a fixed density c — 
Eq. (|2"5|) suggests that (T)m decreases as ~ 1/M, whereas 
in reality the decrease must be slower for small nucleation 
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FIG. 8: Mean nucleation time vs density at a = 1, b = 3, 
and M = 10 6 . The critical density is indicated by a vertical 
dashed line. 




FIG. 9: (Color online). Largest cluster size versus time for 
a = 0.5, 6 = 1, ioi = 5, c = 0.66, M = 10 5 . Results from 
three independent Monte Carlo runs are shown. 

times (large M) since V(t) — » as t — > 0. In particular, 
the mean-field dynamics suggests that for a wide range 
of M values (T)m quasi-saturates at (T)m ~ Uag, since 
the critical cluster size is almost never reached before 

t tlag- 

VII. SIMULATION RESULTS 

Numerical simulations of the zero-range model show 
clear evidence for the existence of a metastable state prior 
to condensation. In Fig.Elwe show the largest cluster size 
as a function of time for three separate Monte Carlo runs 
in the case a = 0.5, b = 1, w\ = 5, M = 10 5 . For each 
run the system was started in a random uniform initial 
condition with density c = 0.66 (for these parameters 
c cr ~ 0.56). It is clearly seen that after a short equilibri- 
ation period the system fluctuates in a metastable state 



T 




20 40 60 80 100 



n 

FIG. 10: Distribution over small cluster sizes in metastable 
state: results for Monte Carlo simulation (crosses) compared 
to prediction of J^J with (w) = w nar and n cr calculated nu- 
merically from 1151 and 1171 (dashed line). 



before a condensate appears. The critical cluster size is 
observed to be around 400 in good agreement with the 
prediction n cr ~ 330 from Eqs. (I5fr and (fTTjl (see Fig.EJ. 
However, the metastable time is about an order of mag- 
nitude larger than predicted. 

In Fig.EHwe show the distribution of cluster sizes (for 
small clusters) averaged over the metastable state of one 
such run. The distribution is in good agreement with 
Eq. JSJ) with (w) = w ncr , thus supporting the assumption 
of quasi-stationarity. 

In the analytical treatment of the previous section we 
calculated the mean time for the maximum cluster size to 
exceed n cr under the assumption that the current in the 
metastable state is constant. In practice, of course, the 
metastable current also fluctuates (and the fluctations 
are greater when w n depends more strongly on cluster 
size, e. g., for the a = 1 case compared to, say, a = 0.5). 
Simulations suggest that these fluctuations can destroy 
the metastable state in cases where the metastable cur- 
rent w naT is close to the current of the condensed phase 

Woo- 

In contrast, for parameters where the metastable state 
is well-separated from the condensed state we find rel- 
atively good quantitive agreement between theory and 
simulation. For example, in Figs. ITD and ll2l we compare 
the average simulation values of critical cluster size and 
nucleation time to the theoretical predictions for a range 
of densities in the case a = 0.5 and 6 = 3. In these sim- 
ulations we crudely identified the end of the metastable 
state as the point when the current out of the largest 
cluster had been less than the average system current for 
50 consecutive Monte Carlo time steps. 

We find that our mean-field theory fairly accurately 
reproduces the critical cluster size but systematically un- 
derestimates the nucleation time. This discrepancy may 
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FIG. 11: Critical cluster size versus density for a = 0.5, b = 3, 
Wi = 5, M = 10 s (c cr ~ 0.27). Crosses show simulation 
data (averaged over 10 Monte Carlo histories), dashed line is 
prediction of Eqs. 1151 and 11711 . 




FIG. 12: Nucleation time versus density for a — 0.5, 6 = 3, 
w\ = 5, M = 10 s (c cr ~ 0.27). Crosses show simulation 
data (averaged over 10 Monte Carlo histories), dashed line is 
prediction of Eqs. <|2*31 and <|2H . 





FIG. 13: The fundamental (flux-density) diagram for two 
different sets of control parameters: a = 0.5, b = 1, wi = 5 
(top); a = 0.5, 6 = 3, wi = 5 (bottom). The branches of 
metastable homogeneous state are shown by dotted lines, the 
critical densities c cr are indicated by vertical dashed lines. 



correlations, etc., our simulations show that the simple 
mean-field approach provides a good qualitative descrip- 
tion of the metastable state and its dependence on den- 
sity. It thus represents an important first step towards 
more refined theories. 



be partly due to the presence of (weak) dynamical cor- 
relations between the numbers of particles in the boxes 
in the fluctuating metastable state. Namely, the appear- 
ance of a large cluster with n ~ n cr is likely to be accom- 
panied by a slight depletion of the surrounding medium. 
Furthermore, we only calculated the mean first passage 
time and ignored the probability that a cluster reaches 
n cr + 1 and is immediately driven by a fluctuation back 
below n cr . Monte Carlo histories which involve such a 
fluctuation back into the metastable state before a con- 
densate is established, would increase the average simu- 
lation nucleation time above the theoretical prediction. 
Despite the neglect of current fluctuations, dynamical 



VIII. CONCLUSIONS 

In this paper we have studied a zero-range model mo- 
tivated by the observation of phase separation in traffic 
flow. In particular, we have demonstrated the existence 
of a metastable state which exists prior to condensation 
and obtained analytical estimates for the critical cluster 
size and mean nucleation time. This new insight into the 
mechanism of cluster formation in the ZRP is supported 
by the simulation results discussed in the previous sec- 
tion. In the context of traffic flow we now conclude by 
further exploring the significance of this metastable state. 

The relation between density c and flux j of cars is 
known as the fundamental diagram of traffic flow. The 
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average stationary flux can be calculated as follows 



3 = 2jQ(") w n , 

71=1 



(25) 



where Q(n) is the probability that there is a car in a 
given cell (in the original model) which can move for- 
wards with the rate w n . Note that only those cars con- 
tribute to the flux, which are the first in some cluster. 

oc 

Hence, Q(n) — (pP(n)/ P(m), where ip is the frac- 

tion of cells, which contain such cars. This fraction can 
be calculated easily as the number of clusters divided 
by the total number of cells. These quantities fluctuate 
in our model. For large systems, however, they can be 
replaced by the mean values. The mean number of clus- 
ters is equal to the mean number of non-empty boxes 

oo 

M Y^, P( n ) m the zero-range model, whereas the mean 

n=l 

number of cells, i. e., the mean length of the road is 
(L) = M + (N) = M (1 + (n)) = M/(l - c), as we have 
already discussed in Sec. IIIII Hence, Q(n) = (1 — c) P(n) 
and lj2"5"jl reduces to 



3 = 0-- c) (w) . 



(26) 



The mean stationary transition rate (uu) depends on the 
car density c. For undercritical densities c < c cr , this 
quantity is the solution of ijSJl. For overcritical densities 
we have (w) = in the phase coexistence regime, as 
discussed in Sec. II VI therefore in this case the fundamen- 
tal diagram reduces to a straight line 



3 = ~ c) 



C>C a 



(27) 



In the metastable homogeneous state at c > c cr the 
mean transition rate (w) together with the critical cluster 
size n cr can be found from the system of equations (fTSf) 
and (|17fl . which allows calculation of the metastable 
branch of flux j. 

The resulting fundamental diagrams for a — 0.5, 
w\ = 5 and two values of parameter b are shown in 
Fig. E| As we see, the shape of the fundamental dia- 
gram, as well as the critical density and location of the 
metastable branch depend remarkably on the value of b. 
These features will also depend on the values of a and w± . 
The metastable branch ends abruptly at certain density 
above which Eqs. (|15f) and (|17fl have no real solution. It 
corresponds to a relatively small, but finite value of the 
critical cluster S1ZG TliQy . 

Note that a metastable branch is also observed in simu- 
lations of cellular automata with slow-to-start rules 0, 
0,^^. In our examples, however, the metastable branch 
is located at larger densities and decreases with increas- 
ing of c over a certain wide range of values depending on 
b, a and w\ . The simulations of the previous section sug- 
gest that when this metastable branch is well separated 
from the condensed section of the fundamental diagram, 
our picture is robust even in the presence of fluctuations. 



In summary therefore, we believe that by suitable vari- 
ation of parameters our simple model can reproduce some 
important features of real traffic flow. There is much 
scope for further investigation, both analytical and com- 
putational. 



APPENDIX: CRITICAL CLUSTER 
DISTRIBUTION FOR < a < 1 

Here we focus on the cluster size distribution at c = c cr 
in the case < a < 1, presenting a detailed calculation 
of the asymptotic behaviour for large n. 

According to JSJ), at c = c cr where (w) = Woo, we have 



In 



P(0)w c 



- J In (l + drn + C + S(n) , 



where C is a constant, and S(n) — > at n — > oo. The 

latter follows from the fact that each term with m = k 

in the sum generates terms oc k~ y with y > when 

the logarithm is expanded in Taylor series, and for each 

k 

of these terms we have k~ y = J mT v dm + O (fc _2/_1 ) 

k-i 

at k — > oo. It ensures that the difference between the 
integral and the sum in lA.lfl is finite and tends to some 
constant C at n — > oo for a > 0. The remainder term 
5(n) is irrelevant for the leading asymptotic behaviour of 
P(n). By expanding the logarithm and integrating term 
by term, for < a < 1 and a ^ f > g> j> • • > we obtain 



P(n) oc J^J cxp 



fc=i 



{~b) k n 



l-ko 



k 1 — ka 



(A.2) 



where [1/cr] denotes the integer part of 1 /cr. The cases 
where cr is an inverse integer are special, since a term 
oc 1/m appears in the expansion of the logarithm, giving 
rise to a power-like correction to the stretched exponen- 
tial behaviour, viz. 



P(n) 



n exp 



fc=l 



(~b) k n 



l-ko 



1 - ka 



(A.3) 



for cr 



The known result for cr = 1 can 



ill 

2 ' 3 ' 4 ' ■ 

also be obtained by this method: it corresponds to the 
power-like prefactor in (|A.3(1 . Only the linear expansion 
term of the logarithm is relevant at 1/2 < cr < 1, so we 
find P(n) oc exp [— bn 1 ^ a /(l — cr)] in the limit of large 
n. The first two terms are relevant for 1/3 < cr < 1/2, 
the third one becomes important for 1/4 < a < 1/3, and 
so on. Eqs. i|A.2|) and i|A.3|) represent an exact analytical 
result at n — > oo which we have also verified numeri- 
cally at different values of a and b. In this form, where 
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the proportionality coefficient is not specified, Eqs. I|A.2|) 
and (|A.3|) are universal, i. e., they do not depend on the 
choice of w% . 
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